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ABSTRACT 

Evidence from pulsar wind nebula symmetry axes and radio polarization ob- 
servations suggests that pulsar motions correlate with the spin directions. We 
assemble this evidence for young isolated pulsars and show how it can be used to 
quantitatively constrain birth kick scenarios. We illustrate by computing several 
plausible, but idealized, models where the momentum thrust is proportional to 
the neutrino cooling luminosity of the proto-neutron star. Our kick simulations 
include the effects of pulsar acceleration and spin-up and our maximum likelihood 
comparison with the data constrains the model parameters. The fit to the pul- 
sar spin and velocity measurements suggests that: i) the anisotropic momentum 
required amounts to ~ 10% of the neutrino flux, ii) while a pre-kick spin of the 
star is required, the preferred magnitude is small 10 — 20 rad s~^, so that for the 
best-fit models Hi) the bulk of the spin is kick-induced with Q ~ 120 rad s~^ and 
iv) the models suggest that the anisotropy emerges on a timescale r ~ 1 — 3s. 

Subject headings: pulsars: general — stars: rotation - supernovae 



1. Introduction 



Th e high velocity of pu lsars, at least an order of magnitude larger than their parent pop- 
ulation flHobbs et al.ll2005f). h a s been a long standing puzzle. Of the many models proposed, 
all except the classic iBlaauwl (1196 ll ) binary break-up mechanism rely directly or indirectly 
on the large release of gravitational binding energy ass ociated with the core collapse event 



Since the fastest obser ved pulsars have v > 10^ km s ^ (ICordes. Romani. &: Lundgrenlll993 
Chatterjee et al.ll2005l ). i.e. bulk velocities > 3x 10^"^ c, binary break-up is clearly inadequate 



and a viable mechanism must tap an appreciable fraction of the co llapse momentum. How- 



ever, a wide range of physical processes can still be invoked (e.g. iLai. Chernoff. fc Cordes 
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200ll : iJanka et al.ll2005l ). so gross energetic constraints alone do not solve the pulsar kick 
problem. 

Rec ent progress has come f rom studies that probe the vector nature of the kick. For 
example iFryer fc Kusenkd (120061 ) argue that neutrino-driven kicks (which may even be crucial 
for successful explosion) will produce fast supernova remnant (SNR) velocities along the 
pulsar velocity vector, while ejecta-driven (recoil) kicks should have fast ejecta opposite to 
the pulsar motion. The pulsar itself also contains kinematic evidence of the kick geometry, 
since both the velocity and spin directions are in principle measurable. These linear and 
angular momenta are fossils of the core collapse which can remain inviolate for as much as 
10® yr before acceleration in the Galactic potential perturbs the velocity or torques perturb 
the spin. T hus mechanisms that are spin -dependent can leave a fossil record in the kick 



distribution (ILai et al. 



2001 



distribution ( Spruit fc Phinney 19981 ). 



Romanil 120051 ). Conversely, kicks can leave a record in the spin 



The key to making such tests is accurate measurement of pulsar motions and pulsar 
spins. For the former, the most important progress has been the dramatic improvern ent in 
the quality and quantity of VLBI astrometric proper motions (IChatterjee et al.ll2005l ). This 
work makes it possible to have a vector proper motion (and in a number of cases a parallax, 
and hence a 2-D space velocity) for most moderately bright, > few mJy, pulsars within 2- 
3 kpc. An even larger sample of appro ximate proper mot ions has been obtained from refined 
analysis of long term timing records (iHobbs et al.ll2005l ). Together these methods provide 
useful space velocities for several dozen young pulsars within 3 kpc. Association with parent 
supernova remnants can provide a few additional proper motion estimates for interesting 
pulsars at larger distances. 

More difficult to measure is the spin direction. The classic approach is to adopt the 
'rotating vector' (RV) model for pulsar polarization and assign the projected position angle 
(PA) of the spin axis to the absolute (Faraday 'rotation measure' RM corrected) position 
angle of the polarization at the phase of maximum PA sweep. This method suffers major 
systematic uncertainties due to poor matches to the expected RV model sweep as well as 
an intrinsic 7r/2 ambiguit y according to whether the emission is in the norma l or orthogo- 
nal mode; early analyses (jPeshpande. Ramachandran. fc Radhakrishnarull999l ) were not at 
all conclusive. However, recent improved polarization measurements and RM corrections 
have provided convincing statistical evidence for a projected 2-D correlation between the di- 
rection s of polarization PA and the proper motion (jjohnston et al.ll2005l : I Wang. Lai, fc Han 
2006al ). Even better, in several cases excellent CXO images have made it possible to measure 
accurately the pulsar spin position angle and its inclination from t he plane of the sky using 
the symmetry axis of the surrounding pulsar wind nebula (PWN) (INg &: Romanil |2004J ). Fi- 
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nally, in a few cases, one can also independently estimate the original post-supernova spin 
of the pulsar. With the 3-D spin orientation and amplitude, comparison with the projected 
2-D motion allows even more powerful tests of the nature of the birth kick. 

In this paper we explore the possibility of comparing the new vector data on pulsar spins 
and speeds with the dynamics of core collapse. The structure of the paper is as follows. In 
section §2.11 we describe three simple models of the luminosity and size evolution of a proto- 
neutron star, which we will follow in estimating the effects of a birth kick. In section ^2.21 we 



summarize the effect of momentum deposition on kick and spin (following ISpruit fc Phinney 



19981 ) and show how the ii itegrated kick controls the birth properties of the pulsar. We also 



comment on the results of IWang. Lai, fc Hani (j2006bl ). whose related but somewhat simpler 



study of kick evolution appeared while we were writing up this paper. In Section ^ we 
discuss the sample of nearby and/or young pulsars that we will use for comparison with the 
models, emphasizing the important objects for which we have more complete knowledge of 
the post-kick state of the star. Section §1] describes a population simulation that allows us to 
compare the distribution of birth properties with the observed neutron stars for the several 
models. Section ^ summarizes the fit results that constrain the basic parameters in the kick 
models. We conclude by discussing the status of the spin-kick comparison, the insight that 
the present data give us into the physics of the pulsar kick, and the prospects for further 
refinement of these tests. 



2. Core-Collapse Luminosity-Driven Kicks 



In this study, we consider some simple models for a kick associated with the large lumi- 
nosity of the post-bounce proto-neutron star. Our assumptions are: i) that the asymmetric 
momentum (kick) scales roughly as the emergent neutrino luminosity, possibly with a de- 
layed onset; ii) that the net thrust originates from a region fixed with respect to the star; 
and Hi) that the thrust may be at a finite angle to the local normal. Of course, this net 
thrust may simply be the surface average of a complex distribution of local forces. However, 
we assume that the sum of these has a stable component with respect to the star on the 
timescale of t he initial rotation and co oling. If the thrusts are impulsi ve and rand o m then , 
as argued by ISpruit fc Phinneyl (119981 ) and illustrated numerically by IWang et al.l (l2006bl ) , 
the final velocity a nd sp in tend to become randomly oriented. Note also, as discussed by 
Spruit fc Phinneyl (119981 ). that if there is no initial rotation, the result of a single kick is an 
orthogonal spin; we will show that the preferred models have a small initial rotation, but 
the bulk of the spin is introduced by the kick itself. 



The simplest realization of this picture would be asymmetric neutrino emission mediated 
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by strong 5— fields. lArras fc Lai Jl999al jbl) suggest that due to weak interaction parity non- 
conservation, the neutrino emission depends on the magnetic field, both through simple field- 
dependent opacities and through helicity-sensitive scattering in the presence of a strong field. 
They suggest that the total asymmetry in the neutrino fiux is ~ Cli^isE""^ + 0.002i?i5/T 
where B = lO^^Sis G is the magnetic field strength, E^, is the mean neutrino energy in 
MeV and T is the temperature in MeV. Clearly, very large fields are needed to produce a 
significant net kick. However, this model does provide a natural 'anchor point' on the star 
and, through the field angle at the surface, a plausible mechanism for non-normal kicks. 

Other possibilities also exist. For example, the ne t momentum raight be generated by 
asymmetric neutrino coupling to matter in the envelope JSocrates et"aL or by 'neutrino 

star-spots' induced by the magnetic field ( Thompson fc Duncan IQQsh. or even by hydrody- 



namic motions associated with accretion onto the neutron star (e.g. iBlondin fc Mezzacappa 



20061). In the last case, fall-back instabilities may apparently build to large amplitude. 



Scheck et al.l (120061 ) have presented a detailed set of 2-D hydrodynamic simulations suggest- 
ing that large recoil velocities of neutron stars occur naturally when m = 1 instabilities with 
strong narrow down-flows develop. Here the kick velocity is dominated by the gravitational 
attraction of the dense down-flow. (A neutrino hot spot associated with this flow does in 
fact generate anisotropic emission, but its kick, a ~ 10% correction, is opposite that of the 
gravitational term). Our model does not, of course, capture the full dynamical content of 
such hydrodynamical kicks, and certainly the neutrino light curve does not fully describe 
the thrust. However this luminosity proflle does provide a characteristic duration for ener- 
getic processes associated with the core collapse. The physical interpretation is thus not as 
clear as for the direct neutrino kick scenario, but the characteristic energetics and time scale 
modeled in our flts can still probe the underlying physics. 



2.1. Proto-Neutron Star Luminosity Models 

We first start from a simple cooling model, where the neutron star is optically thick to 
neutrinos and the hot envelope is a dominant contribution to the moment of inertia. If the 
star's gravitational binding energy is promptly released and lost as neutrino radiation from 
the surface of the proto-neutron star, then the proto-neutron star is in thermal equilibrium 
and the evolution of its radius R^s 'quasi-stationary' and determined by the conservation 
of energy: 

, , , d GM^ 
Lu{t) 



dt R 



NS 



We sh ow here results using the neutrino luminosity L^, from Fig. 3 in lBurrows. Klein. &: Gandhi 



(1l992h but have tried different Ly{t) from the literature. The results are quite insensitive to 
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the details of any standard cooling curve. In the simulation, the initial proto-neutron star 
radius is set to a generic value of 100 km. As shown in Fig. [1], a final radius of 14.7 km is 
obtained, which is determined by the total neutrino fluence. 

At the opposite extreme, we can assume that both the envelope and core are hi ghly 
degenerate at bounce (for reviews: iJanka et al.ll2004j : iKotake. Sato. &: Takahashill2006l ). In 
this case, R^s is independent of Li^{t) and the kick essentially acts at a constant radius 
throughout the cooling. We call this extreme the 'static' model. In our simulations we take 
this as -RAr5=10 km, constant throughout the kick. 

Results from more physical numerical simulations suggest real neutron star behavior 
is between the two limiting cases above. To illustr a te thi s, we examined the resu lts of 2- 
D hydrodynamic simulations from iRampp &: Jankal ((2002) and iJanka et al.l (120041 ). These 
authors employ a Multi-Dimensional Boltzmann Transport and Hydrodynamics (MDBTH) 
code to follow the neutrino advection and diffusion. In Fig. [2] we reproduce the mass shell 
trajectories and show the i^e photosphere as it evolves during the initial cooling. As one sees, 
the neutron star initially has a large but optically thin envelope. We will assume the thrust is 
imparted near the neutrinosphere, but the whole star is forced to co-rotate. The time varying 
moment of inertia can be followed by integrating across these shells. The simulations using 
this evolution are referred to as 'u -transport' models. Fig. [1] compares the basic behavior for 
the models showing L„{t) and the effective RNs{t)- The 'static' model has, by assumption 
the same neutrino light curve as the 'quasi- stationary' picture, albeit at fixed Rns- In the 
case of the 'u -transport' model, the simulations available stopped at t = 0.46 s (marked). 
We extrapolate with simple matched exponentials to take the curves to 50 s. 



2.2. Kick Dynamics 

We wish to model a (possibly non-normal) net thrust at a fixed point r on the proto- 
neutron star surface which scales with the L^{t) multiplied by a constant asymmetry rj. The 
thrust is given by 

F^t) ~ r^^ . 
c 

The angular momentum J = IQ evolves according to 

J = rxF^dt + liQi , 
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where /j and fij are the initial moment of inertia and angular velocity before the kick re- 
spectively. Similarly, the linear momentum is determined by 

p = j F^dt . 

As we are modeling isolated young neutron stars, three angles determine the kick geometry 
as illustrated in Fig. [31 a is the polar angle of the kick position, 6 is the off-axis angle of the 
kick from the surface normal, and is the azimuth angle of the kick about the normal, with 
= closest to the initial spin axis. Note that we are ignoring any initial space velocity of 
the progenitor, since with typical values of only a few km s~^, it is negligible compared to the 
kick velocity of several hundred km s^^. However, some modest fraction of presently isolated 



young neutron stars are released from the break-up of pre-ex isting binaries (jPewey fc Cordes 



19871 : iTauris fc Takenslll998l : iFryer. Burrows, fc Benzlll998l ). In this case, these may have an 



initial velocity of as much as lOOkms"^ (the Blaauw mechanism). If there is no relation 
between this initial velocity and the initial spin, then this pre-natal velocity may be added 
to the post-kick speed in quadrature. The main effect is to blur the correlations described 
below. If, however, there is a fixed relation between the orbital speed and the pre-kick spin 
or the kick itself, then this fourth angle must be introduced to the geometry. As these 
cases represent a small fraction of the single neutron star population, we ignore them here; 
however such considerations may be important to descriptions of kicks for neutron stars 
presently found in binaries. We defer discussion of such models to a future communication. 

During the kick, directions of r and F^ evolve according to the rotation 

f ' = TZf 
K = T^F, 

where the rotation matrix TZ is given by 

+ (1 — cos (5 — cos 5) — sin (5 Vt^^zi^ — cos 5) — Vty sin 5 

Vtr^Vtyil — cos5) — VtzSmS Vt^ + {1 — Vt^) cos 5 — cos 5) + sin 5 

— cos 5) — fiy sin 5 VtyVt^i^ — cos 5) + Vt^^^^^ VlI + [1 — Vtl) cos 5 

Q, = QQ = Q{Qx,Qy,Qz) is the angular velocity and 6 = QAt is the angle that the star 
rotates in time At. 

Given a set of initial conditions, the above equations can be integrated numerically to 
obtain the kick velocity, initial spin period and alignment angle. Note that in this paper, 
we call the initial pulsar spin period Pq which is the value right after the kick. The pre-kick 
spin is denoted by Pi (fij), which is a specified parameter of the model. 
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Fig. m shows two example integrations with identical initial parameters except for the 
pre-kick angular velocities The kick is against the initial spin with a component along the 
surface, thus the angular velocity decreases. However, with Qi = 30 rad s~^, the right panel 
shows that the torque from the kick is large enough to perturb the spin axis to a new direction, 
where the kick is more effective in spinning up the proto-neutron star. Note that the initially 
fast spinner ends up nearly aligned, while the slow spinner is closer to orthogonal. Thus, 
even though the post-kick spin periods Pq are quite similar in the two models, information 
about the pre-kick spin is encoded in the vector direction of Qq. In practice, the typical 
direction and degree of alignment also depends on the thrust orientation and duration; we 
show below that comparison with multiple observables can separately probe these various 
effects. 



3. Pulsar Sample 



Linear and angular momenta are relics of neutron star birth and, if well preserved, 
can provide important probes of the core collapse physics. Given their vector nature, the 
alignment angle between them is an important observable. As discussed above, this is par- 
ticularly sensitive to the kick timescale and geometry and can thus give unique information 
on the kick dynamics. While objects with measurements of all three observables (velocity v, 
initial spin period Pq and alignment angle 'dn-v) are rare, the correlations between these are 
particularly powerful at testing the models. 

Only recently have improved pulsar observations made such comparison possible. For 
example, VLB I astrometry provides unprecedentedly h igh precision proper motion mea- 
surement (e.g 



flHobbs et al.l 12004 



Brisken et al. 



2002 



Zou et al. 



Chatter jee et al.l 120051 ). Also, pulsar timing observations 



20051 



1) have recently provided additional (albeit often more 



uncertain) proper motion values. Altogether these provide several dozens of useful proper 
motion measurements. For pulsar spins, the direction of the projected spin on the sky plane 
has been inferred from the radio polarization profile by adopting the rotating vector model, 
although uncer t ainty in the emission mode leads to an inherent 7r/2 ambiguity. Recently, 
Johnston et al.l (120051 ) obtain the spin orientations of 25 pulsars from r adio polariza t ion ob - 
servations and noted significant alignment with the pulsar velocities. I Wang et al.l (l2006al ) 
reanalyze archival polarization data and arrive at similar conclusions. Next, the discovery 
of X-ray jets and torus structure in PWNe provides an even more powerful (and largely 
mode l-independent) way t o obtain the pulsar spin orientation for a modest number of ob- 
jects. iNg &: Romanil (|2004l ) developed a method of pulsar wind torus fitting to measure the 
3-D (not projected!) orientation of the spin axis. So for the vector (2-D) speeds and spin 
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orientation, progress has been good. In contrast, pulsar initial spin periods remain very 
difficult to esti mate. Currently, only a handful o f pulsars have initial spin period estimates 
(e.g. Table 7 in iFaucher-Giguere fc Kaspill2006l ). As it happens, most of these have PWN 
structure indicating their spin orientations, but many do not have a precise proper motion 
measurement. In some cases, we can make an estimate of the space velocity from the struc- 
ture of the associated SNR, but we must be cautious in applying such model-dependent 
values. 

To select a sample of pulsars with 2-D velocities, we start from the list of pulsar proper 
motions from iHobbs et al.l (120051 ). selecting young isolated pulsars with > 2a proper mo- 



tion measurements. A few objects are added or update d with new measurements (e.g. 
Chatterjee et al.ll2005l : IZou et al.ll2005l : iNg &: Romanill2006l ). Galactic acceleration can affect 
the trajectories of old pulsars, so we restrict our sample to r < 5 Myr (we also have included 
the interesting pulsar PSR B1133+16, which at r = 5.04 Myr is just above this limit). Also, 
differential Galactic rotation can add appreciably to the observed space velocity of distant 
pulsars. When the precise distance is unknown, this can cause large errors in reducing the 
speed and position angle to the local standard of rest. Accordingly, we restrict our sample 
to (i < 3 kpc for pulsars with velocities from astrometric or timing analysis. Table [T] lists 
the pulsar samples used in this study. The first group of 20 objects have only transverse 
velocity measurements, hence they are called '1-D pulsars'. The next group, referred to 
as '2-D pulsars', have both a velocity and a projected spin alignment angle. For IC 443 
and PSR B1800— 21 the spin orientations are obtained fr om X-ray PWNe, all other 14 spin 



measu rements are from radio polarization observations (jjohnston et al.l l2005l : IWang et al. 



2006al ). As noted, the proper motions, must be corrected to the birth local rest frame for 
comparison with the spin directions, as this may affect the position angle by several degrees. 
Finally, the last group consists of 9 '3-D pulsars' which have estimates in all 3 observables, 
namely the velocities, initial spin periods and alignment angles, thus allowing for a full vector 
comparison. These pulsars place the strongest constraints on the model parameters. 

We have checked the consistency of the three subsets of pulsars, comparing their veloc- 
ity and alignment angle distributions with the Kolmogorov-Smirnov (K-S) test. The K-S 
statistics indicates a probability of 0.34 that the velocities of the '1-D' and '2-D pulsars' 
are drawn from the same distribution. Comparison of the '3-D' samples with the summed 
'1-D' and '2-D' ones gives a similar value of 0.33. For the alignment angle, the probability 
of having the same distribution for the '2-D' and '3-D' pulsars is 0.48. Since the samples are 
relatively small, these tests are not very stringent but at least there is no evidence that the 
different measurement methods introduce significant biases in to the sample properties. 



We next summarize the status of the several of the particularly interesting '3-D pulsars'. 
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Crab Pulsar (PSR B0531 +21) iNg fc Romanil (120061 ) have measured the pulsar's proper 
motion with HST astrome try and compared w ith the spin vector, which was obtained from 
pulsar wind torus fitting (INg fc Romanil |2004| ). The initial s pin period is derived from the 
historical age, measured braking index and spin parameters (IManchester fc Taylorlll977l ). 



Vela Pulsar (PSR B0833—45) Tracing back the motions of dense supernova ejecta 
knots, lAschenbach. Egger. fc Triimpeii (119951 ) inferred a supernova explosion site 14^9 ± 7'.2 
away from the puls ar's current position. With the updated pulsar proper motion from 
Dodson et al.l (120031 ). the kinematic ag e of the system i s then r = 20 ± 10 kyr. Using the 
measured braking index of n = 1.4±0.2 (ILyne et al.lll996l ). an initial s pin period of 30=b20 ms 
is estimated. PA of the pulsar spin axis is obtained by torus fitting (INg fc Romanil 120041 ). 



PSR B 0540— 69 Based on the e stima ted pulsar age from the SNR expansion velocity 
( iReynoldd Il985l ). Manchester et al.l (119931 ) s uggest an initial sp i n per iod of 38.7 ms. The 
velocity of the pulsar has an initial estimate by lSerafimovich et al.l (120041 ) using HST/WFPC2 
images. For the pulsar spin, Ng & Romani (2007, in preparation) fit the X-ray pulsar wind 
torus to obtain its orientation. 



PSR J0538+2817 \Ng et al.l (120071 ) report a VLBA astrometric measurement of the pul- 
sar proper motion. They also discovered symmetric structure around the source in the 
X-ray observation, which, if interpreted as polar jets or pulsar wind torus, indicates the 
spin axis. From th e offse t of the pulsar position with respect to the SNR geometrical cen- 
Romani &: NeJ ( 2003 ) argued for a long initial spin period, which was confirmed by 



ter 



Kramer et al.l (120031 ) with the timing proper rnotion results and strengthened by the precise 
VLBA proper motion measurement (INg et al.ll2007l ). 



PSR B 195 1+32 The pulsar proper motion is determined by VLBA astrometry (IMigliazzo et al 



2OO2I ). Assuming the pulsar was born at t he cent e r of th e shell CTB 80, the authors deduced 
the pulsar age and its initial spin period. iHesterl (I2OOOI ) argue that the Ha lobes observed in 
the HST image correspond to pulsar polar jets, which indicates its spin orientation. Radio 
polarization measurements could provide a valuable test of this hypothesis. 

PSR J0537—6910 The X-ray PWN observation shows bar-like extended structure around 
the pulsar. The symmetric structure, if interpreted as pulsar equatorial outflow, indicates 
the pulsar spin vector. Ng & Romani (2007, in preparation) apply pulsar wind torus fitting 
to measure the PA of the symmetry axi s. The pulsar birth site may be inferred from the 
radio image (Fig. 3 in lWang et al.ll200ll ). where the peak of the PWN radio emission is at 
a = 05^37'^44!5(5) S = -69°10'11"(2) (J2000). Then, adopting the kinematic age of 5 kyr 
( jWang fc Gotthelflll998l ). the pulsar velocity is 630 km s~^. While the initial spin period is 
unknown, with its current P = 16 ms we can safely assume the pulsar had a Pq range of 
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8 — 16 ms. 



PSR J1833-1034 ICamilo et all (120061 ) suggest a young kinematic age < 1000 yr for the 
pulsar, which imphes an initial spin period of Pq > 55 ms. They also estimate the PA of the 
pulsar spin axis from the ellipticity of the PWN. Ng & Romani (2007, i n preparation) f itted a 
model torus to the X-ray data and obtained similar results. Following ICamilo et al.l (120061 ). 
we found the geometrical center of G21.5-0.9 at a = 18^33"33!4(3) , S = -10°34'12"(3) 
(J2000). If this is the explosion site, the offset from the pulsar position gives an estimated 
velocity of 125 km s~^. 

PSR JO205+6449 The CXO images of this pulsar in 3C58 have been noted as having a 
prominent central torus and polar jets, and fitting of these provides a robust spin position 
angle (Ng & Romani 2007, in preparation). The position of the pulsar birth may be inferred 
from the SNR radio expansion center to be at a = 02^05' "31!32(3) , 6 = +64°49^58^7(5) 



(J200 0) (M. F. Bietenholz, private communication). Recently iGotthelf. Helfand. fc Newburgh 
( I2OO7I ) have identified a thermal X-ray shell in 3C58 with a similar center. I f one assumes 



that 3C58 is the remnant of SN 1181, one gets an initial pulsar spin of 60 ms (IMurray et al. 
2OO2I ) and a rather large velocity of ~ 840 km s^^. The X-ray and radio evidence, however, 
suggest an older age of 3000 — 5000 yr. This would imply a shorter initial Pq and a transverse 
velocity of ~ 140 — 230 km s~^. While we plot values for the commonly assumed 825 yr age, 
these values and the pulsar's status as a '3-D object' are certainly open to question. A pulsar 
proper motion would resolve these difficulties. 

PSR B 17 6— 44 A clear equatorial torus and polar jets were found in CXO data by 
Romani et al.l (120051 ). where a robust spin position angle and inclination were fit. This 
pulsar is often assoc iated with the SNR shell G343. 1—2.3; comparing the total energy of 
the SNR and PWN, iRomani et aP J2005h suggested Pq = QO - 80 ms. Also, if born at the 
geometric center of the SNR shell, the pulsar velocity is ~ 650 km s^^. However, the PWN 
morphology and interstellar scintillation velocity limits make such a high speed questionable. 
Uneven SNR expansion could, for example, make the birth site much closer. Again, an 
interferometric pulsar proper motion would be particularly valuable. 

PSR B1800—21 Radio images spanning ~ 10 years give an astrometric proper motion 
measurement of the pulsar with an accuracy of 2.5 mas yr~^ (iBrisken et al.ll2006l ). The 
Chandra X-ray observation shows symmetric PWN structure around the pulsar, which is 
fitted by Ng & Romani (2007, in preparation) to measure the spin orientation. However, no 
initial spin estimate is available for this pulsar. 

PSR J 1124- 59 1 6 While the PWN im age in the initial CXO HRC observation showed 
noticeable ellipticity (IHughes et al.l l2003l ) , this axis was not clear in a subsequent ACIS 
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exposure. A scheduled deep ACIS observation could help to resolve this structure. A PWN 
axis would be quite useful here as a robust proper motion axis and a reasonable age and Pq 
estimate are available from the highly symmetric parent SNR. 



4. SIMULATIONS 

With our pulsar samples defined, we are now ready to constrain the kick models. We 
proceed by making Monte Carlo realizations of a given kick model, drawing pulsar spin 
and kick parameters and kick angles from the model distribution, simulating many pulsar 
realizations (by integrating according to the equations in §2.2p and projecting onto the sky 
plane to produce the observables for comparison with the pulsar sample. The comparison for 
a given set of true pulsars is made using a maximum likelihood estimate of the probability 
of getting that observed set from the Monte Carlo model. 

Of course, any physical model will have a distribution of values for the kick parame- 
ters rather than single values. What we attempt to extract from our comparison with the 
observed pulsars are the characteristic values in these distributions. For the kick geometry, 
we assume no preferential directions in the kick position a and the azimuth angle 0, hence 
cos a and (p are distributed uniformly between -1 to 1 and to 27r respectively. The off-axis 
angle 6 is, however, a more interesting parameter as it is potentially related to the kick 
mechanism. Assuming the kick follows a Gaussian distribution around the su rface normal. 



9 is characterized by the dispersion a in the Gaussian. Many authors (e.g. iHobbs et al. 



20051 ) fit the pulsar velocity distribution with a Maxwellian. We adopt this form for the 
momentum asymmetry rj, with characteristic value r]Q. We also attempted to use a simple 
Gaussian but this did not give a good distribution for the observed velocities; apparently, 
the larger Maxwellian tail is important for matching the observed sample. For the pre-kick 
angular velocity fli, we assume a Gaussian distribution with mean fl. To minimize the fit 
parameters, we take the dispersion as Q/3. The value used fli represents the angular fre- 
quency that the star with the same momentum will have after the proto-neutron star shrinks 
to its final radius. Our toy model approximates the time profile of the momentum kick as 
being proportional to the rapidly falling neutrino cooling light curve. A convenient way of 
introducing a characteristic kick time scale into the problem is by a delay to the 'turn-on' 
of the asymmetry, here parameterized as an exponential with a time constant r. This has 
the added advantage of approximating some reasonable physical scenarios for asymmetry 
growth. However, as the asymmetry occupies a smaller fraction of the momentum flux in 
the light curve tail, the effective t] must grow. To avoid focus on this strong correlation, we 
thus report the equivalent asymmetry f] here, which corresponds to the same total thrust for 
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the case of r = 0, i.e. 

r] J L,{t) dt = r]o J [1 - exp(-t/T)] L,{t) dt . 

In brief, the kick model is characterized by the 4 'fitting parameters': rjo, 0, r and a. 
The thrust is given by 

F,{t) = [1 - exp(-t/r)] T] L,{t)/c 
and the kick parameters are distributed as 



a, ~ isotropic 



dN ^ „2 ^^^f ,„2 /r)„2\ 



2(f7/3)2 

^ OC sin^exp(-^V2a') 



The factor sin 6 above accounts for the area in polar coordinates. In the simulation, the ap- 
proximation sin ^ is taken to simplify the calculations. Also we tabulate the asymmetry 
results as fj, as noted above. 

Given a set of fitting parameters, the kick velocity and spin vectors of A^sim — 2 x 10^ 
pulsars are simulated. As the pulsar radial velocity is not an observable, the data samples 
are essentially only projected 2-D values on the sky plane. Therefore, the simulation results 
have to be projected before comparing to the observations. This is done by specifying 2 
angles: the inclination angle ( of the spin axis to the sky plane, and the azimuth angle 0^, 
of the velocity about the spin vector. While the latter is totally random, C is sometimes 
actually known from pulsar wind torus fitting. In the simulation, each of the simulated 
vectors are projected Aproj — 50 times, according to the measured C (and its error) of an 
individual pulsar in the sample if this is known, or isotropically otherwise. 

In order to compare the model prediction to measurements, we define a likelihood func- 
tion similar to the typical statistic 

4~exp[-(as-eL)V2K)^ . 

For an observable ^ e {v, Pq, ^n-v}, Cbs ^"^^ sample pulsar's measurement with uncer- 

tainty o"^ and ^gjjjj is the j'-th simulation result. For the initial spin periods Pq and velocities 
based on the DM-estimatcd distance, uncertainties in the measurements are dominated by 
systematic errors which do not follow Gaussian properties. In this case, a boxcar function is 
used 

£13 ^ } 1 if ICbs ~ CsimI ^ ^\ 



^ I otherwise 
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For each pulsar in the data sample, the likelihood function is then the product of all its 
observables, averaged over the simulations and projections 

where N = A^sim x ^proj- Finally, the total likelihood function is the product of all the pulsars 
in the data sample 

c=Y[c . 

i 

It is convenient to define a Figure-of-Merit (FoM) function for comparison between different 
models. As usual, it is defined as the negative log of the likelihood function FoM= — log C. 

5. Results 

Figures [SHT] show the projections of the likelihood function of our three kick scenarios 
for the main model parameters. Each panel shows four curves, with the minimum indi- 
cating the best-fit parameter values. The solid curve shows results obtained using all the 
objects in our pulsar sample; this generally has the tightest, best-defined minimum. How- 
ever, as discussed in §31 some of the spin and velocity estimates for certain '3-D objects' 
are quite model-dependent and can be questioned. We therefore also ran the FoM calcula- 
tions with a 'minimal' set of '3-D pulsars': the widely accepted Crab and Vela pulsars and 
PSR B0650— 69, which is simple, but has large statistical error bars. The results are very 
similar to those for the full set, albeit with larger uncertainties. Finally, we ran fits with the 
'1-D' and '2-D' objects together and just the '1-D objects' alone. Again, the fits are broadly 
consistent with the full set's results, but as expected the minima are poorly constrained. 
Thus we confirm that our handful of well measured '3-D objects' dramatically improve our 
knowledge of the birth kick properties, but apparently do not introduce major biases into 
the fits. 

Our likelihood function is unnormalized. Since it is a product over individual objects, it 
is convenient to rescale the values for the curves computed for different sets of pulsars. This 
removes the gross scaling of FoM with number of objects and facilitates cross-comparison. 
In practice the results with different object sets are normalized by dividing all by the cor- 
responding best-fit FoM values for the i/-transport model for the same data set. Note that, 
with the exception of the '1D-I-2D' (dashed line) fit to the static model, all minima occur at 
FoM > 1. Also the '//-transport' model has narrower minima. This means that this model. 
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in addition to being the most physical, provides the best description of the full data set of 
the three scenarios tested. 

It is useful to understand how the '1-D objects' can provide some model constraints, 
i.e. how the pulsar velocity distribution is sensitive not only to the asymmetry 77, but to the 
other kick parameters, as well. To illustrate, consider the case of large a or small fi, i.e. when 
kick-induced spin is important. When the pulsar is spun up, the induced linear and angular 
momenta are always orthogonal. In this case, the rotational averaging is efficient, generating 
many low velocity pulsars. In contrast, for large Vl or nearly normal kicks (i.e. small cr), the 
initial spin sets the averaging axis, and the final velocity is principally determined by the 
kick location. For example, a polar kick results in a fast moving pulsar, while an equatorial 
kick is quickly averaged to give a low final speed. Thus, for fast initial spinners, the kick 
velocity distribution becomes relatively uniform, as our model assumes no preferential kick 
direction. This is unlike the observed velocity distribution. 

Of course, when only velocity values are available, there is a large co-variance between 
the model parameters, leading to a broad minimum in the multi-dimensional (projected) 
FoM curve. In fact, it is the handful of '3-D objects' that break these degeneracies and best 
localize the minima, as the FoM curves illustrate. Even then, there are significant correla- 
tions in the fit parameters, as illustrated by Fig. [HI The points show the loci of the minimum 
in two dimensional cuts through parameter space for the zz-transport model. Similar re- 
sults are obtained for the other models. To be conservative, we quote here the 'projected' 
multi-dimensional errors on all the model parameters. The single-parameter errors are, in 
general, significantly smaller. Table [2] reports the best-fit para meters for the models, wit h 



la confidence intervals estimated from the bootstrap percentiles ( lEfron fc Tibshiranilll993l ). 
We generated 10^ bootstrap samples and quote la errors from the standard normal 68% 
interval of the bootstrap replica distribution. The two dot sizes in Fig. |H] show the loci of 
the 'lex' and '2cr' minima. 

There are some robust trends in the fitting results. All four of our basic model param- 
eters are significantly constrained and, with the exception of r in the 'static' model, they 
have significantly non-zero values. Thus, we conclude that in all of our scenarios a finite 
pre-kick spin and a finite tangential component to the kick are required. In some sense, we 
consider the zz-transport model to be the most physical. The best parameters for this model 
are closer to those of the static model than to the shrinking 'quasi-stationary' model. Thus, 
the kick is somewhat insensitive to the initial envelope shrinkage of the z/-transport model. 
However, it should be noted that the zz-transport model strongly prefers a finite r ~ 3s. This 
ensures that the kick is insensitive to the initial large envelope phase. Interestingly, if r is 
forced to be small for this model the best fit f/, a and VL move much closer to the preferred 
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'quasi-stationary' values (Fig. [8]). In this limit, the kick is dominated by the large-radius 
initial phase of the cooling curve. 

The final angular velocity (after the proto-neutron star cooling and kick) has a very 
broad distribution, which peaks at ~ 120 rad s^^ for the zz-transport model. Similar results 
are obtained for the other models. With the pre-kick spin f2j of only 10 — 3 rad s~^, this indi 



cates that most (> 90%) of the spin is induced by the kick, as suggested by lSpruit &: Phinney 



(119981 ). Unfortunately, so few initial spin estimates are available at present that the post-kick 



spin distribution cannot yet be compared with the data. 

The distribution of final velocities for the best-fit i/-transport model constrained by the 
full data sample is plotted in Fig. [9l The observed velocities represent reasonable draws from 
that model, with a K-S probability of 0.3, although the model predicts some additional low 
velocity pulsars. Similarly, Fig. (TU] shows the model and data distributions of the alignment 
angles. The left panel compares the model angle distribution and the data. This is the 
alignment angle range that is actually used in the fits. The K-S statistic gives a moderately 
poor fit with a probability of only 0.08 that the data and model have the same distributions. 
Note, however that the '2-D pulsars' with angles determined by polarization have, by fiat, 
^ < 45°. In the right panel we fold the model curve back to i} < 45°. The data match 
this distribution very well, with the K-S probability increased to 0.98. If our model is 
correct, about 20% of the '2-D' values should be > 45°. Note also that a few of the 
'3-D' measurements have large errors and values > 45° are allowed. We have tested this 
by randomly drawing ~3 of the '2-D pulsars', fiipping the table values to = 90° — t?, and 
re-fitting. The majority of the draws produce, as expected, best-fit parameters close to our 
adopted values, but deeper and tighter minima. 



6. Discussion 

We have explored the effect of an idealized model of a pulsar kick: a single finite duration 
thrust applied at a random but fixed point on the surface of a rotating, cooling proto-neutron 
star. This is a reasonable direct description of a class of models where neutrino anisotropy 
is imposed by a structure (e.g. magnetic) locked in the underlying neutron star. It is also 
a reasonable stand-in for other physical models where e.g. convection is controlled by a hot 
spot on the star surface. We note that our modeled momentum thrust can be considered 
as the net anisotropy of any number of kicks, as long as the anisotropy maintains a fixed 
pattern and follows the evolving core luminosity. Our models do not di rectly address the 



effects of multiple, completely random thrusts, but as first described bv ISpruit &: Phinney 



(Il998l ) and as confirmed in recent model integrations of IWang et al.l (l2006bl ). multiple kicks 
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tend to wash out any anisotropy; the net effect of a large number of kicks is to produce little 
or no correlation between motion and spin, contrary to recent observations. Our si mulations 



also do not directly test the gravity-driven kicks described by IScheck et al.l (120061 ) since for 
a rotating star these are almost certainly directed at the pole. Their scenario also does 
not include the kick-induced spin and rotational averaging described here. The technique 
described here for confronting parameterized models with the pulsar data set can, and should, 
be applied to formulations that probe such kick physics. However, we expect that the basic 
parameter estimates here, describing kick direction timescale and magnitude, can be a useful 
guide for evaluating a wide range of models. 

While our model is certainly a simplified parameterization of real physics, it does pro- 
duce some surprisingly complex behavior that can be discerned in the real pulsar data. 
The most interesting effects are only visible in the '3-D pulsars', objects for which we have 
multiple observables and can measure their correlation within the population. Consider, for 
example, the projected angle {^q.^ vs. the post-kick initial pulsar spin frequency Qq (Fig. fTTj) . 
The model distribution shows that at intermediate initial spin periods 40 — 60 ms, we expect 
a substantial fraction of pulsars to be misaligned, with -^q.^ > 45°. This is because during 
the kick, the spin was slow enough that the initial impulse induced the bulk of the spin, at 
orthogonal angles. In contrast, short spin period pulsars (Pq < 20 ms, > 200 rad s~^) 
tend toward alignment due to spin averaging, and very long period pulsars {Pq > 100 ms, 
Qq < 60 rad s^^) are also aligned. The latter effect is 'survivorship'; to retain a slow spin 
even with a decent kick the pulsar must have started as a slow spinner and have been kicked 
nearly radially. Similar correlations exist with the pulsar velocity, e.g. the fastest moving 
pulsars are expected to be slow spinners (Fig. [12]): if the initial kick produces a large spin, 
then continued thrust 'averages away' and the pulsar never reaches high speeds. 

While considering the 2-D correlations, it is interesting to look at the alignment- velocity 
correlation (Fig. fT^ . Here we may include many of the RV-model pulsars. It is not clear that 
these objects follow the expected correlations. First, there appears a dearth of misaligned 
pulsars at slow speeds. However, two effects may explain this. First, truly slow movers 
are selected against in the data sample as they are unlikely to have either a significant 
astrometric or timing proper motion. Second, we should remember that all RV pulsars have 
been forced to have -d^.^ < 45°. It is possible that at least a few of the slower pulsars are in 
fact orthogonal. What about the pulsars at large space velocity in this plot? Well, the model 
suggests that these should show good alignment. Two of the '3-D pulsars' PSRs B1706— 44 
and J0205+6449 in 3C58 show surprisingly large "^q.^ for their speed. However, note that 
the speed plotted f or PSR J0205-I-6449 assumes a birth date of CE 1181. If, as appears 



increasingly likely (iGotthelf. Helfand. fc Newburghl 120071 ). the pulsar is actually > 3,000 



years old, then the object falls near the peak of the expected distribution. Similarly, the 
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large velocity inferred for PSR B1706— 44 depends on travel from the present remnant center. 
Scintillation measurements suggest a much slower speed, requiring an asymmetric SNR or 
even non-association; these again put the pulsar near the peak of the model distribution. 
Further note the large error bar on the alignment angle of PSR J1833— 1034 in G21.5— 0.9. 
This is set by the uncertainty in the SNR expansion center position. It would not be very 
surprising if this object, as well as the neutron star CXOU J061705. 3+222127 in IC 443, 
would appear at ^9^.^, > 45°. 

Recall that for the 2-D polarization objects, set to -^n.y < 45° by fiat, we might expect 
that 3-4 are actually at large i), emitting in the opposite polarizat ion mode. Thi s woul d 



not, of course, significantly dilute the general trend noted by e.g. iJohnston et al.l (120051 ). 
that RV model fits correlate with proper motion vectors. It is not clear if we will be able to 
understand the radio emission well enough to remove this 7r/2 ambiguity. Happily the direct 
X-ray PWN fits do not suffer this problem, so more such measurements should expose a few 
poorly aligned objects. Globally, we expect fast pulsars to be well aligned. Slow pulsars 
will have a significant fraction of poorly aligned objects. Note that additional velocity 
components, e.g. orbital velocity for pulsars released from binaries, will enhance this trend. 

Are there any serious challenges to our proposed model? Well, yes. PSR B1508+55 is 
fast and has a relatively large -^q.^ = 23° ± 7°. This is very difficult to produce with the 
sort of kick scenarios that we describe. We would predict that improved spin measurements 
would prefer smaller i) < 15°. If a large misalignment perseveres, some amendment is needed 
to the kick scenario. 

What about the substanti al velocities which may arise naturally due to convection. 



when m = 1 modes dominate (jScheck et al.ll2006l )? These kicks make vector predictions 



differing from those modeled here. For example, these authors suggest that fast spinners 
produce strong m = 1 modes and thus only these attain high space velocities (opposite to 
the pattern of Fig. [12]). Also, they suggest that slow spinners (< 60 rad s"^) should be 
poorly aligned (cf. Fig. [11]). This picture may also produce a bi-modal distribution of pulsar 
speeds; it is presently unclear if such bi-modality exists, but this would be a important test 
of the model. Perhaps the most interesting vector test of this down-flow kick scenario should 
come from binary pulsars, where a rotation-controlled m = 1 mode predicts a kick along the 
binary orbit normal. We have not discussed binary kicks in this paper, deferring to a future 
publication. 

Also, we should at least comment on the physical plausibility of our fit values. The 
effective asymmetry fj is quite large. If this were strictly induced at the neutrino photosphere, 
this would represent a large asymmetry. In the simple magnetic field picture this requires very 
large (although possibly transient) fields of ~ iq^5.5~w q Thus, a momentum asymmetry 
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mediated by coupling to matter would be quite appealing, as the field requirements are lower 
The pre-kick spin periods are, as noted above, generally quite low. This argues for effective 



core- envelope coupling, possibly magnetic, in the pre-collapse progenitor (ISpruit &: Phinney 



19981). What about the modest a ^ 0.2 - 0.4 ^ 10° - 20° angle of the kick to the normal? 
This does seem quite natural in a magnetically-induced asymmetry model, where a dipole 
centered at ~ RnsO'/2 would cause such an angle. Pictures where the kick is due to matter 
fall-back, would tend to accrete matter with a large specific angular momentum. While 
effective at re-spinning the star, this would tend to have a large a > 1 and would, of course, 
suggest a strong trend to ~ 90°, which are not observed. 

Finally, we come to the time delay in the kick asymmetry. Since the initial neutrino flux 
falls off exponentially, the delay r can be thought of as selecting a characteristic timescale for 
the momentum thrust. As already noted, this has a strong covariance with the kick amplitude 
and with the pre-kick spin with small delays requiring smaller kicks and faster spins. For 
example the best-flt 'static' model is consistent with no delay and has Qi ^ 10 — 15 rad s^^. 
The quasi-stationary model prefers a ~ 1 s delay and has Qi ^ 20 — 40 rad s^^. Interestingly 
all the flts to proper motion data alone prefer some substantial delay, implying that long kicks 
and signiflcant rotational averaging are useful in producing the observed velocity distribution. 
Finally, the i/-transport model very strongly excludes the no-delay solution. As noted above, 
this is required to avoid having the bulk of the momentum thrust while the star still has a 
large radius. What about a physical origin for this delay? In a picture where a large B fleld 
introduces the asymmetry, it is attractive to associate this fleld with early dynamo action 
in the proto- neutron star. Such dynamo theory is not well d eveloped, but it is interesting 



that qualitative discussions (e.g. [Thompson fc Duncanlll993l ) suggest a characteristic time 



of ~ 1 — 3 s for a convectively-driven dynamo to drive the fleld to saturation. If rj grows as 
~ 0.1Bi5{MeV/ E^)"^ , then the kick gains in efficiency as the B field grows to its limit and 
as the neutrinosphere temperature drops. Thus, the preferred fit parameters, a delayed kick 
with modest tangential component and a slow-spinning pre-kick star, are quite consistent 
with our qualitative picture of a magnetic field induced asymmetry of a neutrino-driven kick. 
Only the large asymmetry fraction, implying a very large B field, stretches this scenario. 

As usual, more observations will provide improved tests of this picture. In particular 
we would like to obtain more robust X-ray image measurement of pulsars with low space 
velocities. Improved polarization sweep position angles for a larger data set will help as well. 
However, the real key to progress is accurate measurements of proper motion vectors for 
more young pulsars, especially slow objects. This means that we need improved sensitivity 
to reach to higher distances and long baselines. High bandwidth VLBI studies extending 
over a decade are the key. To push to larger distances, we will also want improved distance 
estimates to allow the measurements to be corrected for differential Galactic rotation. The 
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real challenge in collecting more '3-D pulsar' measurements is the difficulty in estimating 
initi al spin periods. While some estimates may be obtained from SNR and PWN modeling 



(e.g. Ivan der Swaluw &: Wijl200ll ). most will require either historical remnants or very regular 
shells with high precision expansion ages. There are not many such pulsars, and measuring 
braking indices and precise proper motion for these will not be easy. However, with improved 
X-ray measurements for robust orientations, improved radio measurements for polarization 
and motions, and improved population modeling, we can hope for a more refined exploitation 
of this new probe of neutron star birth dynamics. 
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Fig. 1. — Proto-neutron star cooling. Evolution of the stellar radius and the neutrino fluence 
for the quasi-stationary model are shown in the thick solid and dashed lines respe ctively. The 



thin solid and dashed lines show equivalent values for the i/-transport model (jJanka et al. 



20031 ). with the effective i? of a uniform density sphere with the same instantaneous moment 



of inertia. The thin dotted line shows the radius of the neutrino-sphere in the zz-transport 
model. The end points of the numerical simulation are marked; the curves are extrapolated 
to it = 50 s. 
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Fig. 2. — Mass sh e ll traj ectories of a proto-neutron star in the i/-transport model, reproduced 
from I Janka et al.l (120031 ) . Radius of the neutrinosphere is plotted by the dotted line. 
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Fig. 3. — The kick geometric parameters: a is the polar angle of the kick position, Q is the 
off-axis angle of the kick from the surface normal, is the azimuth angle of the kick about 
the normal, is the pre-kick angular velocity. 
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Fig. 4. — Two kick integrations with identical parameters {a — cos~^ 0.6, (f) — 80°, 9 — 
0.2 rad), except fij = 100 rad s~^ {Left panel) and Jlj = 30 rad s~^ {right panel). The sohd 
hnes follow the instantaneous angular momentum, plotting the equivalent angular velocities 
with the moment of inertia at its final value. The dashed and dotted lines are the parallel 
and perpendicular components of the space velocity with respect to the instantaneous spin 
direction. 
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Fig. 5. — Projected FoM for individual parameters of the 'quasi-stationary' model. The 
solid lines show the full sample fit; the long dashed lines show the fit with the '1-D, 2-D 
pulsars' and the 'minimal set' of '3-D objects' (i.e. Crab, Vela and B0540— 69); the short 
dashed lines show the fit with only '1-D' and '2-D' pulsars; the dotted lines show the fit with 
only '1-D pulsars'. Smoothing filters have been applied to the curves and the FoM scale has 
been normalized by the minimum 'i^-transport' value for each data set. For the full sample 
curves, the upper and lower 'Icr' intervals provided by the bootstrap analysis (Table [2]) are 
indicated by the tickmarks. Comparison with the AFoM at these confidence intervals and 
with Table [2] allows the reader to judge acceptable parameter ranges. 




6. — Same as Fig. [5], for the 'static' model. 
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Fig. 8. — Correlations between fitting parameters for the best-fit //-transport model with 
full data sample. Large and small dots are simulations within 1 and 2a of the best fit FoM. 
Note that significant correlations remain between the parameters and that our full projected 
errors are (conservatively) larger than the single parameters uncertainties. 
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Fig. 9. — Velocity distributions of the best-fit i/-transport model with the full data set, 
compared to the observed pulsars with proper motion measurements of > 2a significance. 
The '1-D pulsars' are shaded. 
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Fig. 10. — Alignment angles of the pulsar samples compared to the best-fit z/-transport 
model with the full data set. The '3-D objects' and the 'minimal set' are shaded in light 
and dark gray respectively. Right panel: Simulations are folded around 45° to illustrate the 
effects of the 90° ambiguity in alignment angles measured by radio polarization. 
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Fig. 11. — Distribution of alignment angle ^9^.^, vs. post-kick initial spin for the best- 
fit //-transport model with the full data set. The simulations are projected for ( = 60°, a 
typical value for the '3-D pulsars'. For these pulsars, measurements with model-dependent 
assumptions dominating the errors are shown by plotting the corresponding error flags in 
dotted lines. 
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Fig. 12. — Distribution of pulsar post-kick initial spin vs. velocity for the best-fit //-transport 
model with the full data set. The simulations are projected for ^ = 60°. The velocities of the 
'1-D pulsars' are marked at top. Again, systematic-dominated parameters of '3-D objects' 
have dotted-line error bars. 
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Fig. 13. — Distribution of alignment angle vs. kick velocity for the best-fit //-transport model 
with the full data set. Here we can plot many '2-D pulsars', so the simulations are projected 
isotropically. The velocities of the '1-D pulsars' are plotted on the top. '2-D objects' are 
plotted by the thin lines. Measurements of '3-D objects' have thick-lined error bars, with 
dotted lines for those parameters dominated by model-dependent systematics. 
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Table 1. Pulsar sample used in the analysis. 



Pulsar i;(kms-l) Ai;^(kms-l) )?n.„ (°) Po (ms) C(°) Ref. 



1-D Pulsars 



B0114+58 


190 


±63 










riUloO+o7 


300 


±68 










B0355+54 


61 


+12/-9 










B0402+61 


653 


±121 










B0450+55 


175 


+20/-19 










B0458+46 


72 


±33 










B0540+23 


215 


±78 










J0633+1746 


128 


±3 










B0656+14 


60 


+7/-6 










B0736-40 


238 


+200/-23 










B0834+06 


157 


+20/-19 










B1702-19 


330 


±152 








1 


B1749-28 


630 


±280 








1 


B1818-04 


129 


+21/-21 










B1848+13 


204 


+25/-25 










B2020+28 


256 


+114/-61 










B2021+51 


115 


+18/-15 










B2022+50 


244 


±33 










B2217+47 


337 


±74 










B2224+65 


1605 


+193/-188 










-D Pulsars 














B0628-28 


318 


+61/-64 


5±4 






2 


B0740-28 


259 


+190/-149 


7±5 






2 


B0823+26 


189 


+55/-34 


21 ±7 






3 


B0835-41 


170 


±30 


13 ± 11 






3 


B0919+06 


506 


±80 


32 ± 17 






3 


B1133+16 


639 


+38/-35 


22 ± 2 






2 


B 1325-43 


597 


±254 


31 ±22 






3 


B 1426-66 


150 


+40/-24 


5 ±9 






2 


B 1449-64 


219 


+55/- 18 


1 ±3 






2 


B1508+55 


1082 


+103/-90 


23 ±7 






3, 4 


B1642-03 


160 


+34/-32 


26 ±5 






2 


B1800-21 


347 


+48/-57 


7±8 




90 ±2 


6 


B1842-I-14 


512 


+51/-50 


5 ±15 






2 


B1929-(-10 


173 


+4/-5 


16 ±2 






2 


B2045-16 


304 


+39/-38 


3±6 






3 


IC 443 


250 


±50 


45 ± 10 






5 


-D Pulsars 














J0205-I-6449 


838 


±251 


21 ± 10 


60± 10 


91.0 ±0.2 


6 


B0531+21 


140 


±8 


26 ±3 


19 ±2 


61.3 ±0.1 


6 


J0537-6910 


634 


±50 


3±5 


12 ±4 


92.8 ± 0.8 


6 


J0538-I-2817 


407 


-I-116/-74 


12 ±4 


139 ± 2^ 




6 


B0540-69 


1300 


±612 


34 ±33 


38 ±2 


93.7 ±5 


6 


B0833-45 


61 


±2 


10 ±2 


30 ±20 


63.6 ±0.1 


6 


B1706-44 


645 


±194 


35 ± 10 


70± 10 


52 ±2 


6 
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Table 1 — Continued 



Pulsar 


V (km s ^) 


Av"^ (kms-i) 


(°) 


Po (ms) 


C{°) Ref. 


J1833-1034 


125 


±30 


16 ± 15 


58 ±3 


85.4 ±0.3 6 


B1951+32 


273 


±11 


18 ±5 


27 ±6 


~ 90 ± 30 6 



''The errors in velocity are compiled from either the uncertainties in proper motion mea- 
surement, or distance estimates, whatever larger. 

''An uncertainty of 20 ms is used in the calculation, as to represent a group of objects with 
slow spin. 



References. — Unless specified otherw ise, al l pulsar velocities are f romlHobbs et al.l ||2005I) 
and re feren ces therein. ( 1) Zou et al.' ('2005* ): (2) I Johnston et al.l l|2005h : (31 IWang et all 
l l2006al') : (4) IChatteriee et al.. (,2005.1 : (5) iGaensler et al.l l|2006h : (6) -see text. 



Table 2. Best-fit parameters of tlie kick models. Tlie best-fit values are obtained from the 

lightly smoothed curves in Fig. [5] — [71 







V 


n 


T 


a 


Quasi-stat 














Full Sample 


0, 


.05 


+0.01 
-0.02 


26+11 




39+0 " 


Minimal Set 


0, 


,05 


+0.04 
-0.03 






5+°-5 

U.O_Q 2 


1D+2D 


0, 


.13 


+0.03 
-0.10 




4+2 
^-3 


^■"-0.17 


ID 


0, 


,24 


+0.03 
-0.21 


131^ 


c+a 


5+°-^ 


Static Model 














Full Sample 


0, 


.07 


+0.03 
-0.02 




1 0+2-1 

-■^•^-0.9 


'-'•'->^-0.18 


Minimal Set 


0, 


,07 


+0.04 
-0.03 




0-8ta73 


n Qo+O-28 

"••J^-0.14 


1D+2D 


0, 


,22 


+0.04 
-0.12 




1.4l?| 


'J-O-0.17 


ID 


0, 


,26 


+0.20 
-0.17 


7+3 
' -5 


O q+1.9 

Z.0_2.i 


n 4+0-5 

^■^-0.26 


i^-transport 














Full Sample 


0, 


,13 


+0.02 
-0.06 




o q+0.4 
'J-'J-1.5 


'-'•^^-0.03 


Minimal Set 


0, 


,13 


+0.04 
-0.08 






U-^^-0.03 


1D+2D 


0, 


,13 


+0.10 
-0.06 


7+27 
' -3 




'-'■^'-'-0.03 


ID 


0, 


,17 


+0.16 
-0.05 


8tl 




n 1 



